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ABSTRACT 

In this talk, we briefly review the influence functional path-integral treatment of quan- 
tum Brownian motion. We report on a newly derived exact master equation of a quantum 
harmonic oscillator coupled to a general environment at arbitrary temperature. We apply 
it to the problem of loss of quantum coherence. 
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INTRODUCTION 


Recently there has been considerable interest in quantum Brownian motion. It was 
motivated by possible observation of macroscopic effects in quantum systems. Among 
them are quantum tunneling with dissipation [1], loss of quantum coherence due to system- 
environment interaction [2], just name a few. The newest application of quantum Brownian 
motion is in quantum cosmology, where the issue of quantum-to-classical transition of an 
open system is very important [3]. These issues also appear in semiclassical theory of 
early universe in which noise, fluctuation and dissipation play important roles in particle 
production, back reaction, phase transition, inflation and galaxies formation [4]. In these 
problems, the interaction between a system and its environment is quite complicated giving 
rise to nonlocal dissipation and colored noise. 

The effect of nonlocal dissipation and colored noise in quantum Brownian motion is 
an outstanding problem, which has been studied only to a limited extent. In some limiting 
cases, the quantum master equation (the time evolution equation) for the reduced density 
matrix of the Brownian motion has been derived before by different authors with different 
methods. These cases are all in the class of ohmic environment, for which the dissipation 
is always local [5]. It corresponds to having a linear damping force proportional to the 
velocity of the Brownian particle classically. The noise associated with the dissipation is 
colored at low temperature. 

Our contribution reported in this talk is the derivation of an exact master equa- 
tion for the reduced density matrix of a Brownian harmonic oscillator linearly coupled to 
a general environment (with a general thermal bath spectral density) at arbitrary tem- 
perature [6]. In our model, the environment is a set of bath harmonic oscillators with 
different natural frequencies. The environment is at a thermal equilibrium state. The 
system (Brownian particle) is brought to contact with this thermal bath. The derivation is 
done from first principles of statistical and quantum physics with Feynman path-integral 
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method and Feynman-Vernon influence functional formalism [7]. This master equation 
can accommodate all possible forms of the nonlocal dissipation kernel and nonlocal noise 
kernel. It is a linear partial differential equation with time dependent coefficients. The 
non-Markovian character resides in these coefficients. In particular we examine the cases 
of ohmic, subohmic and superohmic environment and compute these time dependent co- 
efficients numerically. We show that all the previous master equations obtained otherwise 
are just special examples of our master equation. 


INFLUENCE FUNCTIONAL 


Let us briefly review the Feynman-Vernon influence functional formalism of quantum 
open system. Consider a Brownian particle with mass M = 1 and natural (bare) frequency 
ft. The environment is modeled by a set of harmonic oscillators with mass m n and natural 
frequency u) n . The Brownian particle is coupled linearly to each bath oscillator with 
strength C n . The total action of the combined system plus environment is 


S[x,q] = S[x] + Ss[q] + •S'j„t[x,g] 

=1 *{r 2_ 5 fiV } + / 


(i) 


+ 


l ds ^2{C n xq n | 


where x and q n are the coordinates of the particle and the n-th bath oscillators. 

It is well known that the time evolution of the total density matrix of the system plus 
environment p{t) is governed by the following quantum Liouville equation 


d 


ihj t p(t) = [H, p{t )] 


( 2 ) 


In the coordinate representation, the solution of the above quantum Liouville equation can 
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be written as 


+ 00 -j-oo +OO +OO 


where 


p(x f ,qf,x' f ,q' f ;t) = j dx { J dx\ J dq t j dq\ 

— OO — OO —OO —oo 

X J(x f ,q f ,x' f ,q' f ,t | p{x { , qp, X ?•; 0 ) 

•7 (x f 5 9 /> *^/) 9 /» ^ I i ) Qt i x ,• , ) 0 ) 

?n/ ^n/ 


(3) 


(4) 


= Jdx J Dx' J Dq J Dq exp — | 5 [x,g] — S[x',g']| 

Xi x \ »»•' 

is the propagator of the total density matrix in path-integral form. Here q represents the 
full set of bath oscillator coordinates and the subscript i and / denote the initial and final 
variables. 

We are only interested in how the dynamics of the system (the Brownian particle) 
under the influence of the environment (all bath oscillators). The quantity containing this 
information is the reduced density matrix of the system 

-f-oo 4 oc 

p r (x,x')= J dq J dq' p(x,q]x' ,q')S(q - q) (5) 


— OO — OO 


which is propagated in time by the the evolution operator 

4oo 4 oo 


p r (x f ,x' fy t) = / dx i j dx\ J r (xf,x'j,t | x;,x',0) p r {xi,x\^ 0 


-CO — OO 


If we assume that at t = 0 the system and the environment are uncorrelated 


( 6 ) 


p(t = 0 ) = p s ( 0 ) x p e ( 0 ), 


(7) 


then 


Xf Xj 

J r (x f ,x' f ,t |xj,x-, 0 ) = J Dx J Dx' exp ^| 5 [x] - 5 [x']| F[x,x'] 


( 8 ) 
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The functional factor F[x,x'\ in (8), called influence functional, is defined as 

-foe +OO +OO Qf < 7 / 


F[x,x'} = J dqj J dqi J dql I D *J D *' 


— oo — oo — oo qi q[ (9) 

x exp — {s&lg] + Sint q] - Sjfa'] - Sin<[x',g']}p{,(gi,$J,0) 

It is first introduced by Feynman and Vernon [ 7 ], 

For the problem described by ( 1 ), the influence functional can be computed exactly. 
The result is: 

t 3 1 

F[x,x'] = expj - ^ J ds\ J ds 2 x(si) - x'(si)j *7(^1 - s 2 ) jx(-s 2 ) + x'(s 2 ) 


0 0 

t 3 1 


( 10 ) 


where 


- I j dsi J ds 2 x(si) - x'(sj)j u(s! - s 2 )|x(s 2 ) - x'(s 2 ) | 
0 0 

+00 

K*) = I du> I(lu) coth ^/3hu> costas 


is the noise kernel, also 


and 


Vi*) = Ts lis) 


+ OO 

7(5) = J du 


I(u) 


COS LOS 


U) 


(ii) 


( 12 ) 


(13) 


(14) 


is dissipation kernel. Here I(lo) is the dissipation spectral density defined as 

i(u) = v«( w - Wn )-^L 

„ 2 m n u n 

Tl 

The kernels 77(5) and 1/(3) are generally non-local. There exists an important relation 
between the noise and dissipation kernels, known as the fluctuation-dissipation relation. 
It can be written as 

+00 

u(s) = J ds' K(s — s') 7 (s') (15) 


357 


where the kernel I\(s ) is 


+ OC 


duo . 1 

— to coin -pnu) cosujs 

7T 2^ 


K(s) = / 

0 

which is independent of the dissipation spectral density I(u>) 


(16) 


EXACT MASTER EQUATION 


The detailed derivation of the exact master equation of a quantum harmonic oscillator 
with influence functional (10) has been published in Ref. [6]. Here we just give that 
equation below 

*4 + -*“)} 

+ hfi 2 (#)(x 2 — a:' 2 ) pr(x,x' ,t) 

£ 

- ar(i)(i - *■(*.*'.«) (it) 

— iF(t)h(t)(x — x') 2 p r (x,x',t) 

d d 


+ hT(t)f(t)(x — X '){^Q X + g x i) Pr(x,X ,t) 


where the time dependent coefficients are 

2u 1 (i) 

^ 2 (t) = d 2 (t) - 2r(<)« 2 (0 

f /.x _ «<*«(<) , fgft) -Ci(<) 

/W_ ua(0) + 2r(<)u a (0) 

, . ei(0 — C2(0 

/i(0 = u 2 (t) f(t ) + 8 ctii(f) H 

The time dependent functions Ci(t ), d,(s) and ei(f) in (18) to (21) are 
t t t 

c,(0 = / ds\ J ds 2 j ds 3 Tf(t - 5i)[Gi 2 (5i,5 2 ) + U 2 i(s 2 ,5i)]i/(s 2 - S 3 KO 33 ) 


(18) 

(19) 

( 20 ) 

( 21 ) 


( 22 ) 
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(23) 


t 

di{t) = 2 J ds r](t — s)ui(s) 


t{ i(t) = J ds u(s)u i{ s) 

0 

The elementary functions Ui(s) satisfy the following boundary value problem 

d 2 Ui(s) 


ds 2 


+ 2 J ds'rj(s — s')ui(s') + WjUi(^) = 0 


ui(0) = 1, «i(<) = 0 1 x 2 ( 0 ) = 0, u 2 (t) = 1 


and 


(24) 


(25) 


Vi(s) = u 2 (t — s) v 2 (s) = u-i(t — s) 


(26) 


The Green function ^ 12 (^ 1 , 62 ) in ( 22 ) is 


Gi2(si , $ 2 ) 


^l(^l)^2(^2)^(-Sl - -S2) - M2(si)uiCs2)fl(s2 ~^l) 
^1(62)^2(52) - «l(52)«2(52) 


(27) 


A similar expression for <^21 (>§i , ^ 2 ) can be written in terms of Vi(s). In all the above 
equations, the index i runs from 1 to 2 . 

Let us take a closer look of this master equation. The first line corresponds to the usual 
unitary Liouvillian evolution, which is independent of the system-bath interaction. The 
second line corresponds to a time- depen dent frequency shift (frequency renormalization). 
The third line contains a dissipative term with a time-dependent dissipative coefficient 
T(<). The last two lines contain two diffusive terms with time-dependent coefficients. All 
of these terms depend on the system-bath coupling. Further, one can see that all these time- 
dependent coefficients vanish at t = 0 , when the initial uncorrelated condition is assumed 
valid. The frequency shift and the dissipation coefficient depend only on the dissipation 
kernel while the diffusion coefficients depend on the noise kernel. From (11), (12) and (13), 
we find all these time dependent coefficients are determined by the dissipation spectral 
density (14). 
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A simpler closed formula for the time-dependent coefficients in the master equation 
can be found in the weak coupling limit (up to the first order in the coupling constant 
between the system and the bath), 


t 


<5Q 2 (t) = 2 J ds r](s) cos Qs 
0 

(28) 

t 

r(i) = — I ds T](s) sintls 

(29) 

* j 
0 

■i 


Z 

T(t)f(t) = f ds 1 y(s) sinfls 

j C J 

(30) 

0 

* 


= J ds u(s) cos Sis 

(31) 


o 


EXAMPLES 

An important class of dissipation spectral density is 

I(u) = \u>(%) n - 1 e- ! £ (32) 

7 T iO 

where A is the physical cutoff frequency and uj is another frequency scale usually taken to 
be A. The environment is classified as ohmic if n = 1 , as supra-ohmic if n > 1 or as sub- 
ohmic if 0 < n < 1 . It is important to introduce the physical cutoff frequency because, on 
physical grounds, one expects the spectral density to go to zero for very high frequencies. 
It is clear to see that after introducing the physical cutoff frequency, the dissipation kernel 
(13) is a non-local kernel even for the ohmic environment. 

We have numerically computed the coefficients of the master equation given by (18) 
to (21) for three different environments, namely, ohmic (n = 1), subohmic (n = 0.5) 
and superohmic (n = 3). The damping constant is 70 = 0.3 and the cutoff frequency 
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is A = 2000. The bare frequency 12 is determined from the renormalized one, namely, 
_ q2 _|_ ^2 j ( see explanation in Ref. [6]). Both high temperature region 

(T = 10 5 ) and low temperature region (T = 10) have been studied. The numerical plots 
of these time dependent coefficients and detail analysis could be found in Ref. [6]. 


APPLICATION: QUANTUM DECOHERENCE 


As a simple application,, we discuss the damping of the interference between two 
Gaussian wave packets [8]. 

Let ^1,2 (x,t) be the wave functions of Gaussian wave packets located initially (f = 0) 
at x = ±xo respectively with the same initial spread a 

^i t2 (*,0) =Ne~^£ L (33) 


Let ip{x, t) be the wave function of a system consisting of the superposition of these two 
wave packets, 


The density matrix of the system can be written as the sum of three parts 

p(x,x',t ) = 1p(x,t)^(x',t) = p 1 (x,x',t ) + p 2 (x,x',t ) + Pint(x,x' ,t) 
The probability density function 

P(x,t) =| %}>{x,t) \ 2 = p(x,x,t ) 
can also be written as the sum of three parts 

Pl(x,t) = Pi(x,t) + P 2 (x,t) + Pint(x,t) 

By using the influence functional (10) and the master equation (17), we get 

Cl, 2(1,0 = W(t) ex p[- ~ X 2^(1)^ 


(34) 

(35) 

(36) 

(37) 

(38) 
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and 


P t ni(x,t) = 2v / Pi(£,*) v / P 2 (:M) e D(f) cos 


(39) 


where 


and 


Xo(t) = 


»i(0) 

«2(0) 


x 0 


a(t) = 


«i(0) 2a n (t) 

O • 1 / r\ \ • 


1 


1/2 


■ 1 / 2 ( 0 ) ' a 2 ul(0) ' 4cr 4 i/|(0). 

are respectively the position and spread of the wave packet at time t, 


(40) 


(41) 


m = 


2x 0 [u 2 (0)x + mi(0)x 0 ] 


1 + 8an(t)<7 2 + 4ii 2 (0)<7 4 
is the oscillatory angle (which is present even in the absence of the environment) and 

4a n (t)xl 


(42) 


D(t) = 


1 + 8an(t)(T 2 + 4t/j(0)er 4 


(43) 


is the decay factor (which is present only because of the environment). It is this last term 
depicting the decay of interference between the two wave packets which is usually regarded 
as providing a measure of decoherence. 

We have numerically computed the decay factor e~ D ^ for all the cases described in 
the previous section. The results could be found in Ref. [6]. 
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